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Abstract. We consider a model of one dimensional spin-1 bosons with repulsive 
density-density interactions and antiferromagnetic exchange. We show that the low 
energy effective field theory is given by a spin-charge separated theory of a Tomonaga- 
Luttinger Hamiltonian and the 0(3) nonlinear sigma model describing collective charge 
and spin excitations respectively. At a particular ratio of the density-density to spin- 
spin interaction the model is integrable, and we use the exact solutions to provide an 
independent derivation of the low energy effective theory. The system is in a superfluid 
phase made of singlet pairs of bosons, and we calculate the long-distance asymptotics 
of certain correlation functions. 



On the spin-liquid phase of one dimensional spin-1 bosons 



2 



1. Introduction 

Spinor Bose gases exhibit a much richer variety of interesting macroscopic quantum 
phenomena than spinless bosons. A mean field analysis of three-dimensional spin- 
1 Bose-condensed gases [1,2] and exact diagonalization of the spin Hamiltonian for 
this system [3] showed that the ground state can be either ferromagnetic or spin- 
singlet, depending on the sign of the coupling constant for the exchange part of the 
interaction. In an optical lattice Mott phases occur, featuring an interesting behaviour 
in the magnetic sector [4, 5]. Experimentally, the ground state of such gases, its 
magnetic properties and low temperature dynamical properties have been studied for 
^^Na, which has a spin-singlet ground state [6], and for ^''Rb where the ground state is 
ferromagnetic [7-9]. The creation of a one- dimensional (ID) regime for quantum gases 
by tightly confining the motion of particles in two directions offered new possibilities 
for studying macroscopic quantum effects, and a number of advanced experiments have 
been done for spinless ID bosons [10, 11]. 

In this paper we consider a ID liquid of spin-1 bosons with a short-range interaction 
containing a repulsive density-density term and a magnetic exchange interaction. 
We demonstrate that for the case of antiferromagnetic exchange the low energy 
effective action decouples into two parts: the standard U(l) Tomonaga-Luttinger action 
describing density fluctuations, and the action of the 0(3) nonlinear sigma model 
describing the spin sector. The latter model has a singlet ground state separated from 
the first excited triplet by a gap. Therefore the system of spin-1 bosons interacting via 
antiferromagnetic exchange constitutes a perfect spin liquid where spin-spin correlations 
decay exponentially at large distances. We also identify the "order parameter field", 
which has correlations exhibiting a power law decay at T=0. This field consists of pairs 
of bosons 



Hence the fiuctuation superfluidity of spin-1 bosons is really a superfluidity of pairs. 

We derive our results by two means. For weak interactions we perform a 
semiclassical analysis of the bosonic path integral. For a special point where the coupling 
constants for the density-density and exchange interaction are equal to each other we 
study the Bethe ansatz solution. The latter result is of a certain interest since it turns 
out that the integrable model of spin-1 bosons may serve as an integrable regulator for 
the 0(3) nonlinear sigma model. 

2. Low-energy effective action. Semiclassical derivation. 

Our starting point is a three-component Bose gas with Hamiltonian density 



A = - 1/2*^ 



(1) 
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Here ^{x) = "^q^x), is a three-component Bose field and S" are spin-1 

matrices. We now introduce a number-phase representation 

^^{x) = ^/^n,{x)e"f"^^^\ (3) 

with a real unit vector field ncr{x) 

E-^ = i- (4) 



In the number-phase representation the Hamiltonian density takes the form 7i = T-Cq + V, 
where 
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In order to minimize the potential we eliminate no via the constraint (jl]) and define 
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The new phase fields are taken to have ranges 

< 0+ < vr , < 0_ < 27r . 
In the new variables the potential takes the form 

22, /o 2 2\/2 

n_^_n_ + [2 — — n_){n_^ — 



l! 



n± = 77-1 ± n_i. 
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Minimizing the potential iQ gives the solutions (for gi > 0) 

/i 



A: p = — = po , 

^0 

B : p = — = Po , 
90 



, n+ = 0. 

0. 



(9) 

(10) 
(11) 



From now on we concentrate on solution A. Expanding the potential around the 
minimum gives 



V^ = f(p-po)^ + ^ 



+ ... 



(12) 



where the dots stand for terms of higher order in powers of n±,(j). Carrying out the 
analogous expansion for the kinetic energy part ([5]) of the Hamiltonian density gives 

^2 1 



Ho 



Po 
2m 



(9.0+)2 + —{OA-f + (g.r^„)2 ^_ 
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^mpo 



-{d.pf 



(13) 
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Putting everything together we arrive at the following expression for the action 
S= j dtdx [i^^dt^ - n] 



dtdx 



-pdt4>+ - pon+n^dt4>^ - y0n„9jn„ + ...-H 



(14) 



We observe that p, n_|_ and correspond to massive degrees of freedom. In the next 
step we integrate these out, which in our quadratic approximation can be done by 
"completing the squares". At energies well below the gaps for excitations in the p 
and sectors, retardation effects can be neglected and the low-energy effective action 
remains local in time. Integrating over (p — po) generates a contribution 

l^(5^0+)^ 

to the action, while integrating out n+ and generate respectively 
1 



(15) 



and 



1 



1 



1 - — 

2 



4^1 

Finally we parametrize 

and arrive at the following form for the low-energy effective action S = Sq + Sniam 



So 



dtdx 



2go 2m 



(16) 
(17) 

(18) 
(19) 



S, 



nlam 



J dtdx{^^{idtey + smH{dt(l).y) 



Po 
2m 



(20) 



The action 5*0 describes a free massless boson whereas Sniam is a parametrization 
for the 0(3) nonlinear sigma model. Indeed, defining a three-dimensional unit vector 
field m? = 1 by 

sin 6 cos (, 

(21) 



m 



sin 9 sin ( 

cos 6 



we obtain 



Sniam = ^ J ^^^^ [{dtmf - vl{d^mf] , 
where the spin velocity is given by 



Po9i 



m 



(22) 



(23) 
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The action in the sector can be brought to a more standard Luttinger hquid form 

(24) 

where the charge velocity v^. and Luttinger parameter are given by 



5*0 = I dtdx 
2n 



Vr. 



.c=J^, K^ = nJ^. (25) 

We note that for equal coupling constants go = gi the spin and charge velocities coincide, 
i.e. Vs = Vc- 

2.1. Integration Measure 

In the above derivation we have disregarded the integration measure in the path integral 
of our bosonic theory. We still have to verify that it correctly produces the integration 
measure for the nonlinear sigma model (l20l) . We start by considering the Jacobian of 
the number-phase parametrization 

d{^l^.) ^ _4,^^^_^^2_ (26) 
d(p,ni,n_i,0<^) 

Changing variables to n± gives a factor 

-(n^ — n\)dnj^dn_ (27) 

in the integration measure. Expanding n+ around the minimum of the potential V 
leaves 

ntdn^ = 2^/2 sin^ 6 cos OdO (28) 
The integration measure then looks: 

V2 dp d<p+ d<i)^ dn+ dO d^ sm^OcosO. (29) 
When integrating out the (f) field we wish to change variables to 



n^\ll- — = <p. (30) 



In terms of the ^-parametrization f[T8|) this gives 

rf0 = , (31) 

V 2 sin ^ cos 

which turns the measure into 

dp dn^ d(f) d(f)^ [sinOdO d(j)_] . (32) 

The 6, 0_ piece is indeed the correct integration measure for the nonlinear sigma model. 

We conclude that at low energies the spinor-Bose Hamiltonian ([2]) is described by a 
spin-charge separated theory of a free boson IHM describing collective superfluid pairing 
fluctuations and the 0(3) nonlinear sigma model (122|) describing spin excitations. 
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3. Bethe ansatz solution 

We now turn to an entirely different analysis of the model ([2]), which will lead to the 
same conclusion. The first quantized form of the Hamiltonian ([2]) is 



1 ^ (92 1 

where S are spin 5 = 1 operators. It was shown in Ref. [12] that the model fl33l) 
is integrable along the line Qi = go- A simple understanding of certain features of 
the exact solution can be gained by considering the case of two bosons in the strong 
coupling limit (in which we neglect the kinetic energy). The potential energy of the 
pair of bosons at the same position is obtained by adding two spins S* = 1 to form total 
angular momentum J: 

Ej = go + gi[h{J + l)-S{S + l)], 5=1, J = 0,1, 2. (34) 

The singlet J = is always lower in energy than the triplet and quintet states and at 
gi > ^70/2 its energy becomes negative. In this regime the formation of J = bound 
states becomes energetically favourable compared to having scattering states of two 
bosons. As a result the ground state is made of paired bosons. The integrable case 
go = gi lies well within the parameter region where this occurs. 

Let us now return to the Bethe ansatz analysis. According to Ref. [12] the Bethe 
ansatz equations are given by 

N M 

e'"'^ = - n e4(A:. - k,) J] e-2{k, - A,) , 

j=l a=l 
N M 
-1 = ]^e_2(Aa - /Cj) ]^e2(Aa - Afe) , 
j=l 6=1 
^ N N 

E =^E^^ ^ = E^- (35) 

i=l i=l 

where N is the total number of particles, M = [N — S^), E and P are respectively the 
energy and momentum and 

x + icn/2 

en (a:) = : — jr- 36) 

X — icn 1 2 

The parameter c is related to the couplings (70,1 by 

c = mgo = mgi. (37) 

In order to proceed, we now assume that the system is at a finite temperature T ^ /i, 
where the chemical potential /i is the largest energy scale in the problem. As was shown 
in Ref. [12], in this case [13] only three types of (string) solutions of the Bethe ansatz 
equations ( !35l) contribute to the thermodynamics of the model 
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• Bound states of /c's and A's {k-A strings [14-16]): 



cf =K± ic/2 + 7f 
L± _ A _u ;^ /o I x± 



Xf = A^± ic/2 + 6f , (38) 

where the string centres Aj are real, 7, 5 oc exp(— const L). The zero temperature 
ground state is formed by a partially filled Fermi sea of these bound states. 

• Real fc's. 

• A-strings [14-17]: 

= x^n) ^ .^^^^ ^ l)/2 - j] + 0(exp(-LAr)), J = 1, ...n , (39) 

where the string centres aI"'' G M and A/" is related to the number of excitations in 
this sector. 

Substituting (j38l) into (j35l) and dropping the exponentially small corrections 7, S 
we obtain the following set of equations 

^2iLA. = -Q ^(A^ - A^.) JJ p(A, - A;,) , (40) 
j=i p=i 

Nr Np Mr 

e'Lki = X[e^{ki - k,)\[V{h - Ar^\{e.2{h - K) , (41) 

j=l p=l a=l 

Nr Mr 

-1 =l[e.2{Xa-kj)l[e2{Xa-Xb) , (42) 



j=i b=i 



E 



1 

2m 



Aft A',_ 

,2 



Nb Nr 

P = ^2Ai + ^fcp. (43) 



j=i 1=1 
Here = A^"^ + 2A'6, = A"^ - and 

£:(x) = e6(x)e4(x)e_2(a;), P(x) = e5(x)e_i(x). 

We now take the logarithm of the Bethe ansatz equations for the k-A strings forming 
the ground state and then express them in terms of a counting function [16, 18] 



■,(A) = 2A + -^«(^)+<)(^)-»(^ 

4eK^^)-K^^)- 

p=i 

Here 6{x) = 2arctan(x). In terms of the counting function the Bethe ansatz equations 
(gOD read 

y(A,) = -jf- , (45) 

where /-,■ are integer (for Nb odd) of half-odd integer (for Nb even) numbers. We note 
that the momentum is expressed in terms of the Ij as 

r, A^t Nr 

P = ill'^ + Y.^v (46) 

i=l p=l 
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As shown in [12] the ground state of the system corresponds to choosing = M^. = 
2Nb = N and 



+ J 



J 



This immediately tells us that the "Fermi momentum" is 

ttN tt 



2L 



2N, 



27r,^ 1, 

= -yK^n, + -J 

We are interested in the limit 

L ^ oo , Nf, ^ oo , — — 

Ij 

In this limit we can turn the sum over j in 
Maclaurin sum formula, resulting in 

'A- A' 



(47) 
(48) 



Po = const. (49) 
into an integral by means of the Euler- 



y{A) = 2A 



dA' 



Nr 



9 



3c 



A- A 



p=i 



5c 



6 2 



2c 
A-k 



A- A' 



<(A') 



c 



(50) 



Here a'^ f^^A) is are the root densities of "particles" and "holes" for k-A strings. They 
are related to the counting function by 
dy{A) 



2n [a;(A) + <(A)]. 



dA L pv 

Taking the derivative of ( |50l) we obtain an integral equation 

Nr 



(51) 



1 1 
a'{A)+a'^{A) = - + {aQ + a4-a2)*(r' +- V a5(A-/cp) -ai(A- A;p).(52) 

TT A L ^ — ' 



p=l 



Here * denotes a convolution 



f*9 



dy f{x - y)g{y). 



and the integral kernels are given by 

, , 1 nc 



2ti x2 + (nc/2)2 
Equation fl5^ can be recast in the form 



(53) 



(54) 



1 1 

(/ + aa) *{I-ai)*a' =-- a'f^{A) + -^a.,{A - kp) - ai(A - A;p).(55) 

A TT L ^ — ' 



p=l 



In order to derive the Bethe ansatz equations fl4Tl) over the ground state of the system 



we need to express Hj^i ^(^j ~ ^p) terms of the root densities for k-A strings. We 



have 



l[V{k-Ap)= exp\^-zY,0{2 

p=i 



p=i 



k-A., 



5c 



exp I —iL / dA(Tp{A) 

\ J — oo 

exp {2TTiLf{k)) . 



9 2 



k-A 



5c 



9 2 



k-A 



(56) 
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Taking derivatives we have 

f'{k) = (ai - as) * cTp = ai * {I - 04) * ap . 

k k 

On the other hand Eq. ( l55l) imphes that 

1 , 1 



a4 , 



p=i 



where G = ai * (/ + 02) ^. The kernel of G is given by 
G(A) ^ 



2c cosh 



ttA 



As a result we have 



k 



2n 2nL 



-T' 



2c 



r * a[. 



(57) 



(58) 



(59) 



(60) 



where the /c-independent piece x is fixed by the requirement that /(oo) = 0, and F is 
an integral operator with the kernel 



r(x) 



arctan ( tanh 



TT 



(61) 



In order to proceed we need to evaluate T * cr^. For the ground state this is simply given 
by 



r * ai 



dk [r(A + A;) -r(A- A;)]a;,(A) , 



where 



a;(A) = - + dA' K{A - A') ^^(A'). 

J-Qn 



(62) 



(63) 



Here we have defined a kernel 

K = aQ + - a2. (64) 

Using the relations ao(A) = cro(~A), limA^oo o"o(A) = ^ and taking into account that 
by virtue of Eq. ( l62l) the quantity E * is an odd function of k we conclude that for 

k 

the ground state configuration we have x = 0- Let us now consider the limit where 

c, \k\ < Qo- (65) 
Eor large x we have 

E(x) — ie""^/^ , x > c, (66) 

which yields 



E * a'u 



l.„K(^)/;.A.-.;,A).A.,„,(^)^ 



(67) 
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For an excited state with Up particles and Uh holes with corresponding spectral 
parameters A^''' added to the ground state distribution of k-A strings as well as A^^ 
additional real /c's, the integral equation for the total root density o"' = + a'f^ reads 

1 

a'{A) = -+ dA' K{A - A') a'(A') 
^ J-Q 

' "p rih Nr 

J2 K{A - Af) - ^ K{A - A^) + 5^(a5 - a,){A - k,) 



1 
L 



(6^ 



.1=1 r=l j=l 

Here the integration boundary Q may differ from the one for the ground state (Qo) by 
a contribution of at most order 0{L~^). The integral equation fl68p can be solved by an 
expansion in inverse powers of L 



a\A) = a',{A) + -a[{A) + 0{L'') 



L 



(69) 



where (Tq is the root density or the ground state (!62|) . Hence the leading contribution 
(in to r * cr|j is given by (1671) . The contribution due to a'^ can be seen to be of 

k 

the form 



1 



const + O ^exp 



c 



c, \k\ < Q. 



(70) 



The constant (i.e. independent of k) contribution is by construction precisely cancelled 
by X ill f jSOl) . so that 



p=i 

+ —0(^6X7 



2c 



(71) 



Substituting this into (j60l) . ( j56l) and dropping the exponentially small (in Q/c) 
contribution leads to 

YlV{k - Ap) = e'^^ e-*^^^^'^^^"'^/^) JJe_4(A; - kp)e2{k - kp). (72) 
p=i p=i 

Using this in the Bethe ansatz equations (14111 and then rescaling variables by 



, c6 ^ cu 
k = — , A = — , 
TT vr 



(73) 



we arrive at the Bethe ansatz equations for the 0(3) nonlinear sigma model [20,21]: 

n 

j=l J a=l 



,iAL sinh C 



6'i - gj + ivr 6'i - Ua - iTT 
^1 - - ^1 - Ua + ivr 



1=1,. ..N„ 



n 



Ua — 6j + ivr 
Ua — 6j — ivr 



Mr 

n 



- ^fe + iTT 

Ua-Uh- ivr 



a = 1, . . . ,M^. 



(74) 
(75) 
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The contribution of spin excitations to the total momentum is given by 

Nr 

Pspm = Aj]sinh^z. (76) 
1=1 

We note that in addition to (1751) we still have "physical" Bethe ansatz equations that 
determine the rapidities of the particle and hole excitations over the sea of k-A strings 
forming the ground state. It follows from the finite-volume quantization conditions (1751) 
that the dressed phase-shifts for scattering of such particles and holes with excitations 
in the spin sector vanish (up to exponentially small corrections in Q/c). Hence the 
quantization conditions in the gapless sector must be of the form 

rip 

e^^^(^p = n Spp(^' - ^i) n ^^^(^' - , J = 1, • • • , , 

k^i m=l 
rip rih 

e^LPiA^) = II S,,{A1 - Al) n Shh{A1 - A^) , I = 1, . . . ,n,. (77) 

k=l mj^l 

Here -P(A) is the dressed momentum for k-A strings and Sab{A) the dressed scattering 
phases for particles and holes. They can be expressed in terms of the solution to the 
integral equation fl68|) following standard methods [16,22]. 

3.1. Dressed Energies 

We will now show that the energy in the spin sector is given by 

Nr 

Espin = At;,^cosh^z . (78) 
1=1 

Our starting point are the equations for the dressed energies for k-A strings e'(A) and 
real k's K{k) [12] 

1 2 J2 rQo 



A^ r° 

e'iA) = — - — - 2/i + / dA' K{A - A') e'(A' 
m Am J^Qn 



Qo 

t.2 rQo 

K{k) = n+ dA{a5-ai){k-A)e'{A). (79) 



Here /i is the chemical potential related to Qo by the requirement 

e'iQo) = 0. (80) 
Using the integral equation for e' in the equation for k, we find that 

POO 

K{k)= dA[G{k- A) + G{k + A)]^{A), (81) 

where G is given in Eq. ( 1591) . In the limit c, |A;| ^ Qo this simplifies to 

K{k) ^ cosh(—] - dA e-'^^/^ e'(A). (82) 
^ c / cJq^ 
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Expanding e' around Qo we have 

de' 



e'(A) 



dA 



Qo 



(A - Qo) + 



Noting that the Fermi velocity in the charge sector is defined as 
de' 



dp 



dA 
dp 
dA 



(83) 



^4) 



where p{A) is the dressed momentum for k-A strings and using that this may be 
expressed as p{A) = y{A) by virtue of Eqs. (146|) and (l45l) . we conclude that 



de^ 
dA 



and thus 



27rcT'(A) 



K(k) ^ v.coshf—) —a'iQo) e-"'3o/c_ 

\ C / 71 



(85) 



(86) 



Evaluating the right-hand side of Eq. (167|) by using an analogous approximation and 
rescaling k as before we conclude that indeed 



K,{6) = Avc cosh 6'. 



(87) 



4. Scaling Dimensions. 

The analysis of the finite-size spectrum is complicated by the fact that the ground state 
is made of strings. A priori this could lead to complications related to deviations of 
the strings from their ideal forms (see, e.g. [23]). However, as the critical sector of our 
theory is a simple compactified boson we do not expect such complications to play a 
role. The dressed charge Z = Z{Qo) is defined in terms of the solution of the integral 
equation [24] 

Z{A) = 1 + dA' K{A - A') Z{A'). 

J-Qo 

We see that we have 

Z = na'{Qo). 

It is easy to see that the dressed charge only depends on the dimensionless ratio 
c 



(89) 



7 



Po 



At small 7 we find that 



(90) 



(91) 



On the other hand, for 7 — ^> oo the dressed charge approaches 1 from below. This 
implies that Z is not a monotonic function of 7, but instead has a local minimum. The 
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finite size spectra of energy and momentum are expressed in terms of tfie dressed cfiarge 
as [25] 



E{ANi„D,N+,N^ 



2TTVr 



L 

2tt 



(AN,) 
4Z2 



Z^D^ + N+ + N_ 



(92) 
(93) 

Here kp is given by Eq. AA'';, is tlie cliange in bound state number compared to the 
ground state and D is an integer. Tlie corresponding spectrum of scaling dimensions is 
given by 



P{ANb, D, iV+, N^) = — [N+-N^ + D{ANk)] + 2fc^D. 

1j 



h{ANb, D, N+, N_ 



DZ 



DZ - 



(AiVfc 
2Z 

{AN, 
2Z 



1 2 



(94) 



Let us look at some specific examples. The lowest charge-neutral excitations with 
momenta 2kF and respectively correspond to AA*";, = N± = 0, D = 1 and 
D = Nb = N_ = 0, N+ = 1 

/i(0,l,0,0) = ^ , 

/i(0,0, 1,0) = 1. (95) 
The lowest excited state with charge 2 corresponds to AN, = 1, N± = D = 



/i(l, 0,0,0) 



8Z2 



(96) 



5. Correlation functions. 



Let us now turn to the large-distance asymptotics of correlation functions. We start by 
considering the gapless sector of the low-energy effective action, which is given by 



dtdx 



(97) 



Given that the original phase fields 0i 2 were 27r-periodic, the field is in fact vr-periodic 
(see Eq. (|H])). As a result, the spectrum of scaling dimensions of primary operators in 
the Gaussian model (p71) is given by 

1 / 2m \ ^ - 1 / 2m \ ^ 

h{m,n) = - i —= + uVk ] , h{m,n) = -[—=- uVk ] ,(98) 



where n and m are integers. We note that (p8ll agrees with (194|) in the appropriate limit, 
where we have \/K = 2Z. The operators corresponding to these scaling dimensions are 
constructed by introducing the field dual to 0+ by [26] 



71" 

[6+{x), 0+(x')] = i-sgn(x - x') 
The boson number is given by 

A^ = 1 j dx dj+{x), 



(99) 



(100) 
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The local operators corresponding to the spectrum of scaling dimensions fl98l) are 

02m,n = e^'""^++'"'+ ; m,n integer. (101) 

These operators carry charge 2m, reflecting the fact that our fundamental objects are 
pairs of bosons and have two-point functions of the fornjjl 

(O2™,n(x)OL,n(0,0))= [^] e-^'>^, (102) 

where qq is a short- distance cutoff determined by the gapped degrees of freedom we have 
integrated out in order to arrive at the action fl22|) . flMl) . We are now in a position to 
write down expressions for operators that do not involve the spin sector. These are the 
pairing field A and the density, which are given by 

A = - ^2/2 , 

: p : = - po ■ (103) 

As A carries charge —2 we conclude that it takes the form 

oo 

A(X) =P0 ^n02,n{x) , (104) 

n=—oo 

where An are numerical coefficients. We will see below that A2k+i = 0. Similarly, as 
the density is neutral it is given by 

oo 

: p := -dj+ + V BnOo,n{x) • (105) 

TT ^ — ' 

n=— oo 

Like for the pairing field, it turns out that by virtue of the statistics of the Bose field 
all odd coefficients vanish, i.e. B2k+i = 0. Let us now turn to the Bose field itself. At 
low energies the semiclassical analysis gives 

4 . (106) 

Since the spin part of the bosonic operator is proportional to the slowly varying unit 
vector field m, we conclude that the triplet excitations of the nonlinear sigma model 
occur at small momenta g ~ 0. To take into account the quantization of total charge we 
should complete Eq. fll06p with a piece containing the field dual to 0+, which amounts 
to the substitution 



v/;^e*^+ ^ y Po + J2 C2„e2-^++2^nupoXei^+ . (107) 

n 

We note that only even exponentials of the dual field are allowed by virtue of the bosonic 
statistics of the fields which requires the Lorentz spin of the operators appearing in 
Eq. (11071) to be integer. At T = we then obtain 

{T^^i{r,x)^^>{0,0)) = 5,,,, (T,m,(r,x) m,(0, 0)) /charge (r, x) , (108) 

I We note that if were a 27r-periodic field, the allowed scaling dimensions would be h2Ti{m,n) = 
i ^-^^ + 2n^/K^ , /i2,r(m,n) = ^ ~ 2nv^^ and correspond to operators C'm,2n- This is the 
case for Haldane's "bosonization of the boson" [26]. 
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/charge (t, x) — Pq 



Po 



al 









Q2^2mpox 



VT + ix 



h.c. 



VT — IX 



(109) 



The correlation function of the sigma model is well approximated by the expression [27] 

(T,m,(r, x) m,(0, 0)) ^ Zi^o(A + a;Vt;2), (110) 

where Z is a normalization. We note that finite temperature correlation functions can 
be calculated by using the conformal mapping in the gapless sector and the method of 
Ref. [28] for the nonlinear sigma model. 

The expression for the Bose-fields feed back to the low-energy projections of the 
density and pairing field, which are bilinears in the Bose fields. Given that only even 
bosonic exponentials of the dual field appear in the expression for "^cr, the same must 
hold for : p : and A. This requirement sets the coefficients of 02n,2fc+i in the expressions 
for : p : and A to zero. The two-point functions of both these operators undergo a power 
law decay at large distances 

(A(x)At(O)) 



x^ 



+ 2|y42|^cos(27rpox) 



x^ 



Kc+l/Kc 



+ 



(: p{x) : : p(0) :) 



+ 2i?2 cos(27rpoa^) 



x^ 



Kc 



+ . . 



:iiii 



(112) 



2tt'^x'^ 

We see that the oscillating piece of the two-point function of the pairing field A always 
decays more rapidly with distance than the oscillating piece of the charge density. 



5.1. Optical Lattice 

The model ([2]) arises as the low-density continuum limit of a three-component Bose- 
Hubbard system 

Hbh = a] ^a^+i,^ + h.c. + Uo^n'^j + Sj - p^ uj , (113) 

j j 

where Uj = Yla^] a^j,^ spin-1 operators on site j. An interesting 

question is what the phase diagram of the model flllSp looks like. A Density Matrix 
Renormalization Group analysis [4] has established the regions in parameter space where 
Mott phases occur. In order to address this issue by analytical methods one would need 
to bosonize the Hamiltonian (11131) for strong interactions and large densities. Assuming 
that the operator content suggested by Eqs. fll06p and (I107P remains unchanged, this 
would occur when the oscillating part of the Bose field becomes sufficiently relevant so 
that the oscillating parts of the operator turn into a relevant perturbation in the 
entire charge sector. The perturbation is of the form 

5S = \ [ dtdx Oo,iit,x). (114) 
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It is relevant in the RG sense if Kc < ^. It would be interesting to check whether the 
phase boundaries established in [4] coincide with the loci in parameter space where Kc 
becomes 1/2. 

6. Discussion. 

In this work we have analyzed a continuum model of a spinor Bose gas in one spatial 
dimension. Using standard semiclassical methods we showed that for weak interactions 
the low energy degrees of freedom are described by a spin-charge separated theory of a 
free boson and the 0(3) nonlinear sigma model. We then turned to the integrable line in 
the model. From the exact Bethe ansatz solution we obtained an independent derivation 
of the low-energy effective theory. An interesting by-product of our analysis is the 
demonstration that the integrable model provides a simple, integrable regularization of 
the 0(3) nonlinear sigma model. Finally, we determined the long-distance asymptotics 
of certain correlation functions (including the "order-parameter" correlator) in the 
framework of the low-energy description. The dominant fluctuations in the theory 
are of superfiuid singlet pair type. Our analysis for dynamical correlations pertains 
mostly to the weak-coupling regime. The opposite limit of (infinitely) strong repulsive 
density-density interaction should amenable to a treatment along the lines set out in [29]. 
Other interesting open questions are how the ID spinor Bose gas behaves under a 
quantum quench [30] and what happens in the unbalanced case where the densities of 
the three species of bosons are fixed at different values. Finally, it would be interesting 
to generalize our analysis to spin S bosons [31]. 
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